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Off-body pressure, forces, and moments for the Gulfstream Low Boom Model are com- 
puted with a Reynolds Averaged Navier Stokes solver coupled with the Spalart-Allmaras 
(SA) turbulence model. This is the first application of viscous output-based adaptation to 
reduce estimated discretization errors in off-body pressure for a wing body configuration. 

The output adaptation approach is compared to an a priori grid adaptation technique 
designed to resolve the signature on the centerline by stretching and aligning the grid to 
the freestream Mach angle. The output-based approach produced good predictions of cen- 
terline and off-centerline measurements. Eddy viscosity predicted by the SA turbulence 
model increased significantly with grid adaptation. Computed lift as a function of drag 
compares well with wind tunnel measurements for positive lift, but predicted lift, drag, 
and pitching moment as a function of angle of attack has significant differences from the 
measured data. The sensitivity of longitudinal forces and moment to grid refinement is 
much smaller than the differences between the computed and measured data. 

I. Introduction 

A detailed review of the history and state-of-the-art of sonic boom modeling is provided by Plotkin 1 
and literature surveys are provided by Ozcer 2 and Park. 3 Advancements in Computational Fluid Dynamics 
(CFD) tools and their application have resulted in significant improvements to near-field sonic boom signa- 
ture prediction. Methods have been developed to align initial grids with the Mach cone before beginning 
analysis. 4 Hybrid methods have been developed 5- ' where near-body unstructured grid solutions are inter- 
polated to shock-aligned structured grid methods to increase accuracy. Solution adaptive grid methods 8,9 
have been employed to improve the accuracy of unstructured grid methods for long propagation distances. 
Output-adaptive approaches have been applied to sonic boom prediction in two-dimensions (2D) with a dis- 
continuous Galerkin method, 10, 11 a three-dimensional (3D) Cartesian grid method, 12 and 3D unstructured 
grid methods. 13-15 

When CFD validation studies have been performed on complex geometries, the existing methods perform 
much better for the forward portion of the signature than the aft portion. 4, 14 A notable exception is Signature 
47 of the Shaped Sonic Boom Demonstrator (SSBD) program. 16 There are many sources of uncertainty in 
simulation and measurements of off-body pressure signatures. Wind tunnel test sections can exhibit pressure 
fluctuations on the same order as the signature. 1 ' -19 Some of the uncertainty in the prediction of the aft 
signature may be due to uncertainties in the geometry. 3 

Often, near-field sonic boom predictions are performed with Euler methods, so the effect of a viscous 
boundary layer is absent. Output adaptive Euler and turbulent Reynolds Averaged Navier-Stokes (RANS) 
approaches were applied to predict extremely near-field pressure signatures of a plume in supersonic flow. 20 
Differences in signatures between Euler and RANS simulations motivated the use of RANS coupled with a 
turbulence model for the prediction of off-body pressure signatures of a more complex wing body configura- 
tion. 

A low boom configuration must also have good aerodynamic performance to be a practical design. Im- 
provements in aerodynamic performance can result in reductions in fuel required to perform a mission. A 
reduction in fuel usage also reduces empty vehicle weight, vehicle volume, and environmental impact. Re- 
ducing the weight and volume of the aircraft directly aids in sonic boom mitigation. 21 Therefore, accurate 
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predictions of lift, drag, and pitching moment are important to the design of a practical aircraft and are 
included in this study. 


II. FUN3D Flow Solver 

This study utilizes the Fully Unstructured Navier-Stokes Three-Dimensional (FUN3D) simulation frame- 
work. FUN3D is a suite of tools for the analysis and design of a wide range of aerospace vehicles. The 
user manual, an extensive list of references, and example applications are available at the FUN3D websited 
The FUN3D framework contains multiple algorithms to allow for the evaluation and comparison of different 
numerical methods. The following FUN3D description is limited to the methods utilized in this study. A 
finite- volume method 22 is applied to tetrahedral grids. The flow variables are stored at the vertexes or 
nodes of the mesh. At interfaces delimiting neighboring median dual control volumes, the inviscid fluxes are 
computed using Roe 23 or van Leer 24 approximate Riemann solvers based on the values on either side of the 
interface. For second-order accuracy, interface values are obtained by extrapolation of the control volume 
centroidal values, based on gradients computed at the mesh vertexes using an unweighted least-squares tech- 
nique. 22 The reconstruction limiter is based on the framework of Barth and Jespersen, 25 but the min and 
max functions are replaced with a differentiable function. 26 This limiter, denoted hvanalbada, is frozen 2 ' 
after 1000 iterations. 

The full viscous fluxes are discretized using a finite-volume formulation in which the required velocity 
gradients on the median dual faces are computed using the Green-Gauss theorem. The solution at each 
time-step is updated with a backward Euler time-integration scheme. At each time step, the linear system 
of equations is approximately solved with a point-implicit procedure. 28 Local time-step scaling is employed 
to accelerate convergence to steady state. The Spalart-Allmaras (SA) model 29 is used to model turbulence. 

III. Output-Based Adaptation 

Output-based adaptation requires the solution of discrete adjoint equations. 3, 30,31 These linear adjoint 
equations are solved with a time-marching method that produces the same asymptotic convergence rate 
as the flow equations. 28 To improve adjoint iterative convergence, a generalized conjugate residual (GCR) 
scheme 32 is used with the standard point-implicit procedure as search directions. The linearization of the 
turbulence model and eddy viscosity is included in the adjoint residual. The reconstruction limiter is frozen, 
which removes its dependency from the adjoint linearization. 

Venditti 33 described an output-based error estimation and adaptation scheme which utilizes an embedded 
grid. This has been implemented in 3D 34 within a parallel 35 framework. The embedded grid required for 
the Venditti error estimator can be prohibitively large (eight times the original grid) for 3D problems, which 
effectively limits the problem size of the output adaptive method. A single-grid error estimation procedure 
has been described 3 that produces similar adaptation guidance as the embedded-grid procedure for Euler 3 
and turbulent 36 flows. 

Venditti 33 provided a procedure to calculate a new grid spacing request from the adaptive indicator and 
a user specified error tolerance. In this study, this error tolerance is set to half the total remaining error 
estimate for each adaptation. The anisotropy of mesh elements is based on a scaled Mach Hessian, where 
the element size in the smallest spacing direction is dictated by the output-based spacing request. 33 This 
scaled Hessian is expressed as a metric, which is utilized by the adaptive grid mechanics. The selection of 
the Mach Hessian to determine the anisotropy of the metric is arbitrary, but has produced reasonable results 
for problems with shocks and boundary layers. 20 

Parallel grid mechanics described by Park and Darmofal 3,3 ' are used to adapt the grid. The local 
elemental operators of node insertion, node movement, element swap, and element collapse are utilized to 
iteratively drive the edge lengths to unity in the specified metric. The grid adaptation mechanics suffer 
from robustness problems when performing highly anisotropic adaptation near curved boundaries. Also, 
the discretization in FUN3D is more accurate for semi-structured right-angle elements in the boundary 
layer. 38 To address both of these issues, the near-wall boundary layer grid is frozen during the adaptation 
process. This approach is not applicable to problems with initially under-resolved geometry or boundary 
layers because that portion of the grid is held constant during the simulation. However, this approach is a 
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powerful technique for problems with important off-body anisotropic features, i.e., sonic boom propagation. 
Details of the adaptation technique with constant near-wall boundary layer grid are available in Andren and 
Park 39 and Park and Carlson. 20 

The cost function for the adjoint error estimation procedure is the integral of quadratic pressure deviation 
over a surface s in the domain, 



( P- Poo 
V Poo 


2 

els, 


(1) 


where A s is the area of the integration surface. The adaptation is focused on improving the calculation of 
pressure near this surface also referred to as the sensor. Integrating the square of pressure deviation has 
been shown to produce more accurate signatures with less control volumes than integrating the pressure 
deviation directly. 12 In this study, a cylindrical surface, aligned to the £-axis, is defined for s to resolve the 
signature on and off the centerline. 14 The extent of the integration surface can be optionally restricted to 
the interior of a bounding box to focus on a portion of the cylinder surface. 


IV. Model and Experiment Description 

Two cooperative wind tunnel tests were performed by the National Aeronautics and Space Administra- 
tion and the Gulfstream Aerospace Corporation. The model is designed to have a Quiet Spike™ telescoping 
forward fuselage extension 40 integrated into a low boom and low drag design. 41 Separate boom and per- 
formance wind tunnel tests were performed, each with a different model support sting. The boom model 
is smaller than the performance model. This allows off-body pressure measurements at further relative dis- 
tances from the model and performance measurements at a higher Reynolds number. The model support 
stings are tailored to minimize their interference on the data measured in each test. 


V. Off-Body Pressure Prediction 

The off-body pressure test 41 was performed in Test Section 1 of the NASA Langley Unitary Plan Wind 
Tunnel. 42 Test section cross section is nominally 4- by 4-feet. The wind tunnel employs sliding block nozzle, 
which imparts flow angularity in the test section. 42 Techniques for grit application to control boundary layer 
transition and guidelines for flow angularity are provided by Wassum and Hyman. 43 Grit was only used for 
one of the signatures presented here. The model was supported with a curved blade sting 44 that attached to 
the upper surface of the fuselage, Fig. 1. The wind tunnel test and simulation is performed with 3.5 million 
Reynolds number per foot or 0.613 million Reynolds number for the mean aerodynamic chord. 

The van Leer flux function is selected for the boom propagation simulation because the Roe flux was 
not suitable due to the strong expansion around the blade sting. The simulations were run at Mach 1.6 and 
0.25 degrees angle of attack. The wind tunnel walls are not modeled in the simulation; freestream boundary 
conditions are applied to the outer boundaries. 



Figure 1. Boom model with blade sting. 


The initial grid was generated with the VGRID and SSGRID tools. 4,45 The boundary layer grid is 
constant for all nodes within 0.002377 mean aerodynamic chords of the no slip boundary conditions. This 
corresponds to the first 12 nodes of the advancing layer mesh. The initial grid has 3.5 million nodes. 
Adaptation is performed eight times and the final adapted grid has 11.5 million nodes. The concatenated 
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convergence history of the flow solver is shown in Fig. 2 for the initial and adapted grids. The residuals 
of the five mean flow equations are plotted with solid lines and the residual of the SA turbulence model is 
plotted with a dashed line. The flow solver used 2,000 iterations on each grid. The reconstruction limiter is 
frozen after 1,000 iterations resulting in a rapid drop of the mean flow residuals. The adaptation mechanics 
interpolate the solution, which is used for each flow solver and adjoint solver restart. The mean flow 
equations converge approximately ten orders of magnitude on each grid and the turbulence model converges 
approximately nine orders of magnitude. A loosely coupled time advancement scheme is used for the SA 
model because it allows for guarantees on the positivity of v to improve robustness during initial transients. 29 
The adjoint solver utilizes 10 search directions for GCR convergence acceleration. The adjoint equations 
converge approximately eight orders of magnitude on each grid in 500 iterations (not shown) . 



The symmetry planes of the initial and adapted grids are shown in Fig. 3. The intersection of the model 
and the symmetry plane is shown in white. The curved blade sting is clearly seen above the fuselage and 
vertical tail. The initial grid is shown in the left subfigures and the adapted grid is shown in the right 
subfigures of Fig. 3. The normalized delta pressure dp = (p — Poo)/Poo is shown in the upper subfigures. The 
intersection of the cylindrical boom integration surface and the symmetry plane is shown as the thick black 
line in Fig. 3b and Fig. 3d. The sensor cylinder s has a radius of h/l = 2.0, where l = 13.2 inches is the 
reference body length. Only the negative 2 and negative y quadrant is used. The x extent of the cylinder 
is x/l = [2,4]. The forward limit is ahead of the bow shock. The aft limit of the integration surface was 
set to exclude the strong model support pressure disturbances that were not measured in the wind tunnel 
test. The output adaptation technique increased the grid resolution above the model to reduce the estimated 
error in the pressure below and to the side of the configuration, see Fig. 3b and Fig. 3d. The adapted grid 
upstream of the bow shock and downstream of the model is coarser. The initial grid has more high frequency 
content in the middle of the signature on the centerline. This difference is also observed in a close up of the 
symmetry plane near the model, Fig. 4. 

The increase in grid resolution above the model observed in Fig. 3 is also observed in Fig. 4. The 
difference in the high frequency content of the middle of the signal observed in Fig. 3 is more evident in 
Fig. 4. The aspect ratio of the tetrahedra has been reduced and the grid is less structured in this middle 
region of the signature. Grid density has been increased in the wake of the configuration. The pressure aft 
of the configuration is smoother and the expansion at the tail has a different shape. 

A closeup of the forward fuselage of the model is shown in Fig. 5. The model patch boundaries, surface 
grid, and surface contours are also included in Fig. 5. The delta pressure color contour range is adjusted 
to highlight small pressure variations behind the nose boom shocks. The nose boom shocks and expansions 
are targeted by the output adaptive approach above and below the model. The lower nose boom expansions 
are different between the initial and adapted grid. This may be due to the difference in grid orientation. 
The initial grid is stretched along the freestream Mach angle, but adapted grid stretching is determined by 
the local Mach Hessian. The current Mach Hessian recovery technique successively applies a least squares 
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(b) Adapted grid delta pressure contours. 


I 


(a) Initial grid delta pressure contours. 




(c) Initial grid. 


(d) Adapted grid. 


Figure 3. Symmetry plane of the entire near-field boom propagation region. 
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(a) Initial grid delta pressure contours. 


(b) Adapted grid delta pressure contours. 




(c) Initial grid. 


(d) Adapted grid. 


Figure 4. 


Detailed view of the near-body propagation region symmetry plane. 
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(a) Initial grid delta pressure contours. 


(b) Adapted grid delta pressure contours. 



(c) Initial grid. 



(d) Adapted grid. 




(e) Initial grid turbulence model parameter i>. 


(f) Adapted grid turbulence model parameter i>. 


Figure 5. Detailed view of aft fuselage symmetry plane. 


7 of 22 


American Institute of Aeronautics and Astronautics 





gradient technique. 22 An adaptive grid procedure that utilized a Green formulation for Hessian recovery 
produced better resolution of small scale details than an approach utilizing the successive least squares 
approach. 8 

The SA turbulence model variable v is shown in the lower subfigures of Fig. 5. The region where the 
turbulence model is active above and below the model is well resolved by the adapted grid. The boundary 
layer is treated as fully turbulent. The SA turbulence model rapidly approaches the freestream value at the 
edge of the boundary layer. This edge region is targeted by output error estimation techniques. 46 The peak 
levels of the turbulence are higher for the adapted grid. 

The surface pressures are uneven in both the initial and adapted grids indicating that the surface may 
need to be refined. The source of the high frequency disturbance in the signature propagated by initial grid 
is difficult to determine from the surface grid pressures. The pressures are nonsmooth at the edge of the 
boundary layer for both grids. It appears that the change in topology between the unsheared and sheared 
regions of the initial grid and the edge of the highly refined adapted boundary layer is at the same location 
as this nonsmoothness. 

To illustrate the differences in the initial and adapted grid off centerline, a volume slice at x/l = 0.6 is 
shown in Fig. 6. The initial grid is shown on the left half of each subfigure and the adapted grid is shown 
on the right. The intersection of the volume grid tetrahedral elements and the plane is shown in Fig. 6(a). 
Pressure interpolated to the plane is shown in Fig. 6(b). Silhouettes of the wing, tail, and blade sting are 
seen in the center of each subfigure. The adapted grid has higher resolution in near the body and the nose 
boom shocks on the centerline and off the centerline. The grid has been coarsened outside of the nose shock. 



(a) Intersection of volume grid and plane. (b) Delta pressure contours. 


Figure 6. Volume slice at x/l — 0.6. 


V.A. Pressure Signatures 

The wind tunnel runs are listed in the order presented, Table 1. The distance between the pressure mea- 
surement and the x-axis is h/l. The angle between the symmetry plane and the off-track measurement 
location is </>. The measurements are made with probe 2. 41 Off track comparisons will be provided first to 
illustrate the most significant differences between the initial and adapted grids. Run 197 is a repeat of Run 
148 with grit added to force boundary layer transition. All other Runs did not use grit and were allowed 
to transition naturally. All the simulations assumed a fully turbulent boundary layer without transition. 
The predicted signatures are shifted aft by the distance denoted in the following figure legends to align the 
predicted expansion in the middle of the signature to the measured location. A shift is commonly employed 
for wind tunnel comparisons because of slight variations in wind tunnel test section Mach number and other 
types of flow nonuniformity. 
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Table 1. NASA Langley UPWT Run 
descriptions in order presented. 


Run 

h/l 

(/> (deg) 

Figure 

172 

1.201 

26 

Fig. 

7a 

173 

1.230 

38.5 

Fig. 

7b 

174 

1.338 

56.5 

Fig. 

8a 

175 

1.478 

81.6 

Fig. 

8b 

156 

1.730 

25.5 

Fig. 

9a 

157 

1.756 

35.0 

Fig. 

9b 

158 

1.829 

53.4 

Fig. 

10a 

159 

2.002 

83.9 

Fig. 

10b 

140 

0.49 

0 

Fig. 

11a 

133 

1.2 

0 

Fig. 

lib 

148 

1.7 

0 

Fig. 

12 

197* 

1.7 

0 

Fig. 

12 

t Run 

197 has 

grit applied to 

ensure 


boundary layer transition. All other 
cases are free transition. 


Off centerline data for h/l = 1.201 and h/l = 1.230 is shown in Fig. 7. The series of nose boom shocks 
are better resolved by the adapted grid. Comparisons are made even further from the aircraft centerline in 
Fig. 8. The nose boom shocks predicted by the adapted grid are absent in the initial grid. The largest delta 
pressure is well predicted by simulations on both grids. The adapted grid predicts a similar lowest delta 
pressure to the initial grid near the centerline. Away from the centerline, the adapted grid has a lower delta 
pressure which compared well with the measured signature. 

Off track delta pressure comparisons for h/l = 1.730 to h/l = 2.002 are shown in Fig. 9 and Fig. 10. 
The required shift increases with off track angle. A constant value of 0.0035 was added to wind tunnel 
measurements for Run 159, Fig. 10b, so that the delta pressure would fluctuate around a zero mean upstream 
of the boom shock. Reference pressure can drift before and during measurements. 4 ' Again at this location, 
the adapted grid predicts a lower delta pressure than the initial grid for the expansion, especially for the 
most off track location. 

The predicted and measured signatures for centerline at h/l = 0.49 and h/l = 1.2 are shown in Fig. 11. 
Due to limitations in the model traversing mechanism at this closest location, only the aft portion of the 
signature was measured in Run 140, Fig. 11a. The initial grid produces high frequency signal content for 
x/l = [0.8, 1.0] in Fig. 11a and x/l = [1.7, 1.9] in Fig. lib. The propagation of these peaks can be seen in 
the Fig. 3a, Fig. 4a, and Fig. 5a. The high frequencies of the aft portion of the adapted signature are also 
damped. The output adaptive method is resolving the signature over the entire lower cylinder quadrant and 
the initial grid is designed to primarily resolve the centerline. 

The effect of tripping the boundary layer in the experiment is shown in Fig. 12. The approximate locations 
of the grit applied to the model 48 is denoted as the thick lines in Fig. 13. The open circle of Fig. 12 are the 
natural transition measurements and the filled squares are the measurements with grit applied. The increase 
in the high frequency signal at x/l = 2.4 in Fig. 12b coincides with the approximate location of the grit strip 
on the fuselage. The highest pressure measurement at x/l = 2.6 is closer to the predicted signatures. The 
low pressure region at x/l = [2.8,3] has also been affected by the change in the boundary layer state. The 
second low pressure valley at x/l = 2.9 has been increased as a result of grit application, which improves 
its correlation with the predicted signature. The measured peak at x/l = 3.1 has been reduced, which also 
improves the correlation with the predicted signature. 
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Initial, shift=0.30 

Adapted, shift=0.30 
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(a) Run 158, h/l = 1.829, <j> = 53.4. 



Initial, shift=0.33 

Adapted, shift=0.33 — 
UPWT, dp=0.0035 


2 2.5 3 


x/l 

(b) Run 159, h/l = 2.002, <j> = 83.9. 

Figure 10. Pressure signature far off centerline. 
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Figure 12. Pressure signature on the centerline at h/l = 1.7, with (Run 19T) and without (Run 148) grit to trip the 
boundary layer in the wind tunnel. 
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VI. Force and Moment Predictions 


Force and moment tests 49 were performed at the NASA Glenn 8- by 6-foot Wind Tunnel. 50 The simula- 
tion conditions are 1.56 Mach and 4.82 million Reynolds number per foot or 3.38 million Reynolds number for 
the mean aerodynamic chord. The hvanalbada limiter is used with the Roe flux function and SA turbulence 
model. The drag integration surfaces are shown in gray in Fig. 14. The red model support sting is omitted 
from the simulation force calculation. The wind tunnel walls are not modeled in the simulation; freestream 
boundary conditions are applied to the outer boundaries. 



Figure 14. Gray force integration surface and red model support sting. 


Three different grids are used for this simulation. They contain 2.85, 5.40, and 7.67 million nodes. The 
patches and sources had to be modified to produce the two finest grids. 45 Producing a family of uniformly 
refined unstructured grids remains an open topic. To provide an example of the typical iterative convergence 
for these grids, the iterative convergence at 3.917 degrees angle of attack is shown for all three grids in 
Fig. 15. The continuity equation is plotted with solid lines and the SA turbulence residual is plotted with 
dashed lines. The limiter was frozen at iteration 1,000 and a total of 3,000 iterations were performed. The 
flow equation residual was reduced six orders of magnitude and the turbulence model residual was reduced 
five orders of magnitude from its peak. The convergence rate of all three grids is very similar. 



Figure 15. Convergence history at a. = 3.917 for three grid resolutions. 


The longitudinal forces and moments of the model (excluding the sting and base) are shown in Fig. 16 
and Fig. 17 as a function of angle of attack. Lift as a function of drag and pitching moment is shown in 
Fig. 18. All three grid levels are shown. The simulation was performed at the same angles of attack as 
the wind tunnel measurements. The slope of the predicted lift curve is higher than the measured lift curve 
slope, Fig. 16a. The predicted drag values are larger than the measurements for the entire angle of attack 
range, Fig. 16b. The slope of the predicted pitching moment is steeper than the measurements, Fig. 17. The 
predicted and measured pitching moments agree best near zero lift. The measured change in slope near five 
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(b) Coefficient of drag. 

Figure 16. Longitudinal forces as a function of angle of attack. 
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Figure 17. Coefficient of pitching moment as a function of angle of attack. 


18 of 22 


American Institute of Aeronautics and Astronautics 




Glenn 8x6 



Glenn 8x6 



-0.05 -0.04 -0.03 -0.02 -0.01 0 C 


CM 


(b) Coefficient of lift as a function of pitching moment. 

Figure 18. Lift as a function of drag and pitching moment. 
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degrees angle of attack is also seen in the simulation. The variation with grid resolution is greatest for drag 
at negative angles of attack and large magnitude pitching moment. The variation due to grid resolution is 
much smaller than the difference between the measurements and simulation. 

Lift as a function of drag and pitching moment is shown in Fig. 18. The differences between the measured 
and predicted data in Fig. 16 collapse for positive lift, Fig. 18a. The grid sensitivity is greatest for the negative 
lift portion of the curve. USM3D also showed that the prediction of lift as a function of drag was better 
than lift or drag as a function of angle of attack. 45 USM3D 51 utilizes a cell-centered discretization, so it has 
approximately six times as many control volumes as node-centered FUN3D for the same grid. 

VII. Conclusions 

Simulations were performed to obtain off-body delta pressures, forces, and moments for the Gulfstream 
Low Boom Model. The off-body delta pressure predictions were compared with wind tunnel measurements 
of a model with a blade sting model support. The predicted forces and moments were compared to wind 
tunnel tests of a model with a conventional rear sting model support system. All of the simulations were 
performed with RANS coupled with the SA turbulence model to include the effects of a turbulent boundary 
layer. Off-body pressure predictions have historically been performed with Euler methods. This was the 
first application of a viscous output adaptive technique to predict off-body pressure of a 3D wing body 
configuration. The adaptation mechanics do not currently permit the modification of the grid in the lowest 
portion of the boundary layer. 

The output adaptive approach was compared to an a priori grid adaptation technique that stretches 
and aligns the grid to the freestream Mach cone. The a priori adaptation method resolved small, high 
frequency pressure disturbances on the aircraft centerline, but was less accurate for off-centerline predictions. 
The output adaptive technique simultaneously targeted centerline and off-centerline portions of the signal 
and provided balanced predictions that were good on and off the centerline. If desired by the engineer, 
the centerline resolution can be increased and off-centerline resolution reduced by the tailoring the sensor 
location to only include the region near the centerline. An improvement to the Hessian recovery method or 
an output-based approach that directly determines the anisotropy may improve the output adaptive results. 

One of the presented wind tunnel measurements used grit to trip the boundary layer. The grit produced 
measurable signature changes at the grit location that were not modeled in simulation. The application of 
grit also changed the aft portion of the measured signature, but it was unclear if the fully turbulent boundary 
layer simulations better matched the tripped or untripped case. The eddy viscosity predicted by the SA 
turbulence model increased dramatically with grid adaptation. 

Computed lift as a function of computed drag compared very favorably with wind tunnel measurements 
for positive lift, but predicted lift and drag as a function of angle of attack had significant differences from the 
measured data. Predicted pitching moment as a function of angle of attack or lift had the same general trends, 
but slightly different slopes. The measured and predicted pitching moments were closest near zero lift. The 
forces and moments of the configuration were computed with a series of refined grids. A weak sensitivity 
to grid refinement was noted, but could not justify the differences in the computed and measured data. 
This suggests that the differences between experimental and computational results may include significant 
modeling errors as well as discretization errors. 
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